Analytic Theory of fractal growth patterns in 2 dimensions 



O 

o 
o 

(N 



o 



c3 

■4— » 

c3 



^3 

o 
o 



> 

o 
m 
o 
o 
o 
-i— > 

c3 



i 

O 

o 



X 



Benny Davidovitch and Itamar Procaccia 
Dept. o/ Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 

(version of February 1, 2008) 



Diffusion Limited Aggregation (DLA) is a 
model of fractal growth that was introduced in 
1981 by Witten and Sander Q. It had since at- 
tained a paradigmatic status due to its simplicity 
and its underlying role for a variety of pattern 
forming processes including dielectric breakdown 
|^], viscous fingering [Q, electro-chemical depo- 
sition dendritic and snowflake growth 0, 
chemical dissolution [Q, geological phenomena [H 
and certain biological phenomena such as bacte- 
rial growth [pi and viscous fingering through gas- 
tric mucin pOQ . The algorithm begins with fix- 
ing one particle at the center of coordinates in 
(i-dimensions, and follows the creation of a clus- 
ter by releasing random walkers from infinity, al- 
lowing them to walk around until they hit any 
particle belonging to the cluster. The funda- 
mental difficulty of such growth processes is that 
their mathematical description calls for solving 
equations with boundary conditions on a com- 
plex, evolving interface. Despite tremendous ef- 
forts this and related difficulties defied all at- 
tempts to develop analytic theory of DLA. In 
fact, even the numerical estimates pi]] of the frac- 
tal dimension D of DLA turned out to converge 
extremely slowly with the number of particles n 
of the cluster, leading even to speculations [ [l2| ] 
that asymptotically the clusters were plane filling 
(i.e. D = 2). In this Letter we offer a theory for 
fractal growth patterns in 2-d, including DLA as 
a particular case. In this theory the fractal di- 
mension of the asymptotic cluster manifests itself 
as a renormalization exponent observable already 
at very early growth stages. Using early stage 
dynamics we compute 1.6896 < D < 1.7135, and 
explain why traditional numerical estimates con- 
verged so slowly. The present theory is equally 
applicable to other fractal growth processes in 2- 
dimensions, and we discuss similar computations 
for such growth models as well. 

For continuous time processes in 2-dimensions the 
above mentioned difficulty was efficiently dealt with in 
the past Jl3|-^5[ by considering the conformal map from 
the exterior of the unit circle in the complex plane to 
the exterior of the (simply connected) growing pattern. 
In this way the "interface" in the mathematical plane 
remains forever simple, and the complexity of the evolv- 
ing interface is delegated to the dynamics of the confor- 
mal map. For discrete particle growth such a language 
was developed recently M-19[ , showing that a large va- 



riety of fractal clusters in two dimensions can be grown 
by iterating conformal maps. In this Letter we employ 
this language to show that the fractal dimension of the 
asymptotically large clusters has a surprising and useful 
role as a renormalization exponent in a rescaling theory 
of these clusters in their early growth phases. This find- 
ing allows us to compute the fractal dimension to desired 
accuracy. 

Once a fractal object is well developed, it is extremely 
difficult to find a conformal map from a smooth region to 
its boundary, simply because the conformal map is terri- 
bly singular on the tips of a fractal shape. The derivative 
of the inverse map is the growth probability for a random 
walker to hit the interface (known as the "harmonic mea- 
sure" ) which has been shown to be a multifractal measure 
p0| characterized by infinitely many exponents pl| , p2[ . 
Accordingly, in the present approach one grows the clus- 
ter by iteratively constructing the conformal map start- 
ing from a smooth initial interface. Consider $( n >(w) 
which conformally maps the exterior of the unit circle 
e 10 in the mathematical u>-plane onto the complement of 
the (simply-connected) cluster of n particles in the phys- 
ical z -plane fll6|-|l9|]. The unit circle is mapped to the 
boundary of the cluster. The map <§>( n '(w) is made from 
compositions of elementary maps (f>\ $, 



(1) 



where the elementary map 4>\. g transforms the unit circle 
to a circle with a "bump" of linear size VA around the 
point w — e l6 ' . An example of a good elementary map 
was proposed in fl6[| , endowed with a parameter a in 



the range < a < 1, determining the shape of the bump. 
We employ a = 2/3 which is consistent with semicircular 
bumps. Accordingly the map (u>) adds on a new 
bump to the image of the unit circle under & n ~ 1 \w). 
The bumps in the z-plane simulate the accreted particles 
in the physical space formulation of the growth process. 
Since we want to have fixed size bumps in the physical 
space, say of fixed area Ao, we choose in the nth step 



A« — 



Ao 



|$(n-l)'( e ie„)|S 
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The recursive dynamics can be represented as iterations 
of the map 0a„,0„ (w), 



& n >(w) = <j)\ u Bi ° </>A 2 ,e 2 ° ■ ■ ■ ° 0A»,0„M 



(3) 



The difference between various growth models will 
manifest itself in the different itineraries {9\ ■ ■ ■ 9 n }. To 
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grow a DLA we have to choose random positions 8 n . This 
way we accrete fixed size bumps in the physical plane ac- 
cording to the harmonic measure (which is transformed 
into a uniform measure by the analytic inverse of $(»)). 
The DLA cluster is fully determined by the stochastic 
itinerary {0k}£ =1 . In Fig. |l] we present a typical DLA 
cluster grown by this method to size n =100 000. 




FIG. 1. A DLA cluster, n = 100000. 




FIG. 2. A deterministic cluster with a Golden Mean 
itinerary, n = 100000. 

Other fractal clusters can be obtained by choosing a 
non-random itinerary [ fL9| . A beautiful family of growth 
patterns is obtained from quasi-periodic itineraries, 

9 k+l = 9 k + 2nW , (4) 

where IF is a quadratic irrational number. An exam- 
ple is shown in Fig. |^, in which W is the golden mean 
{VE + l)/2. In Q it was argued that itineraries ob- 
tained from (H) using for W other values of quadratic 



irrationals lead to clusters of different appearance but 
the same dimension, which was estimated numerically to 
be D = 1.86 ± 0.03. In the same paper other determin- 
istic itineraries (not obtained from circle maps) where 
shown to lead to clusters with different dimensions. One 
(trivial) example that is nevertheless useful for our con- 
sideration below is the itinerary 8 k = for all k. Such 
an itinerary grows a 1 dimensional wire of width \/Ao . 

The great advantage of the availability of a conformal 
map is that it affords us analytic power that is not obtain- 
able otherwise. To understand this consider the Laurent 
expansion of (w) : 

= F^w + F (n) + F^w- 1 + F^w- 2 + ... 

(5) 

The recursion equations for the Laurent coefficients of 
$(n)(ty) can be obtained analytically, and in particular 
one shows that fil| 

n 

F[ n) = n[l + A fc ]°. (6) 
fc=i 

The first Laurent coefficient F^ has a distinguished role 
in determining the fractal dimension of the cluster, be- 
ing identical to the Laplace radius which is the radius 
of a charged disk having the same field far away as the 
charged cluster [[17]. Moreover, defining R n as the min- 
imal radius of all circles in z that contain the n-cluster, 
one can prove that Q 

Rn < 4F 1 (ll) . (7) 

Accordingly one expects that for sufficiently large clus- 
ters (to be made precise below) 

F x (n) ~ n^VAo , n->oo, (8) 

as >/Xq remains the only length scale in the problem when 
the radius of the cluster is much larger than the radius 
of the initial smooth interface (which we take as the unit 
circle in this discussion, $(°)(cj) = u>)). 

These observations lead now to the central develop- 
ment of this Letter, and to the most important result. 
Consider a renormalization process in which we fix the 
initial smooth interface, but change Ao, and then rescale 
n such as to get the "same" cluster. Of course we need 
to specify what do we mean by the "same" cluster, and 
a natural requirement is that the electrostatic field on 
coarse scales (i.e. far from the cluster) will remain in- 
variant. In other words, we should require the invariancc 
of the Laplace radius F^ (and possibly of additional 
low order Laurent coefficients) under renormalization. 

Clearly, for a given itinerary {6*fe}^ =1 , is a func- 

tion of n and Ao only. Accordingly, considering Eq.(^|), 
we note that such a renormalization process can reach 
a fixed point if and only if (Ao) attains a nontrivial 
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fixed point function F* of the single "scaling" variable 
x = y/Xon}' . Obviously in the asymptotic limit i>1 
F x (n) (A ) must converge to F* which is linear in x in this 
regime. The main new findings of this Letter are that F* 
exists as a nonlinear function of x, and that (Aq) con- 
verges (within every universality class) to its fixed point 
function F£ already for x <C 1. 

In principle one can demonstrate the convergence to 
F* analytically. This is easy to do in the case of the 
degenerate itinerary growing a wire. In this case we can 
demonstrate convergence after the addition of 2-3 bumps, 
even in the limit Ao — * 0, see Fig. || . For x — ► F*(x) — 
1 ~ x 2 in this case. 
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FIG. 3. The one dimensional wire. We plot Fi(x) — 1 vs. 

1/2 

x — n\ and demonstrate the convergence to the asymptotic 
nonlinear function F*. The values of Ao used are 10 -3 , 10" 4 , 
1GT 6 and 10 -6 . 

For other nontrivial itineraries it becomes increasingly 
cumbersome to demonstrate the convergence by hand. 
With the assistance of the machine we can demonstrate 
the convergence in all the other cases. In Fig. [| we 
present F^ n ' (Ao) — 1 as a function of x for a typical DLA 
itinerary and for values of Ao ranging between 10 -8 to 
10 -3 . We note that for Ao — > the convergence to the 
fixed point function is obtained infinitcsimally close to 
the initial circle for which F^ n ~ ^ =1. In fact, data col- 
lapse (with D chosen right) for this itinerary, as well as 
for all other nontrivial itineraries, is obtained for n > n c 
where n c « 2-k/\/\~q. This is the number of bumps re- 
quired to obtain one-layer coverage of the original circular 
interface. Obviously n\^ D \f\) — > for Ao — » 0, demon- 
strating the convergence to F* for x <C 1. In Fig.|| we 
exhibit the convergence for the Golden Mean itinerary. 
Note that the fixed point functions are different, and they 
both differ from the wire case. The main point of this 
analysis is that convergence to F* can be obtained for x 
arbitrarily small by going to the limit Ao — * 0. 
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FIG. 4. Convergence to F? for the DLA. We plot Fi (x) - 1 
vs. x = n}/ D with D — 1.7. The values of Ao used 
are 10 -3 (upper solid line), 10 -4 (upper dashed), 10 -5 (thin 
solid) 1(T 6 (dotted), 10~ 7 (lower dashed) and 10" 8 (lower 
solid). The fixed point function F* is best approximated by 
the lowest curve in the figure. Convergence to the asymptotic 
nonlinear function F* is seen for smaller x values when Ao 
decreases. 
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FIG. 5. Convergence to F* for the Golden Mean itinerary. 
Data are same as in Fig. 4 but x — n 1 ^ \j/ 2 with D — 1.83. 
The fixed point function F* is best approximated by the low- 
est curve for x > 3 x 10~ 2 where the data with Ao = 1CP 8 
converges to F*. 

The existence of a fixed point function translates im- 
mediately to a calculational scheme. Consider a given 
itinerary {Ok}k=i °f one °f the above classes, and calcu- 
late F x (Ao) for N > n > n c (Ao)- Rescale now Ao — > 

A /s, and calculate F[ n \\ /s) for N > n' > n c (X /s). 
We can compute D from finding the value n' which pre- 
serves the Laplace radius under rescaling of Aq by s: 
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1/D 



(9) 



Since F± is monotonic (as is immediately seen from 
Eq.(0)), there is only one solution to this equation, 



D = 



2[logn' — logn] 
logs 



(10) 



As a first example consider the wire case. Computing 
F^ 10 '(10~ 6 ) we find that for Ao = 10~ 5 the same value 



of F± is obtained for n between 3 and 4. Eq.(10) with 



s = 10 then predicts 0.796 < D < 1.045. Repeating for 



^"^(lO" 6 ) we find the same value of f\ w (10- 5 ) for n 
between 31 and 32. From Eq.@ 0.9897 < D < 1.0173. 

We stress that this precision is obtained when F^ — 1 is 

still 0.0133! Lastly, using i^ (500) (10~ 6 ) = 1.18254 we 
compute 0.9951 < D < 1.0006, and any desired ac- 
curacy can be achieved by increasing n. The reader 
should note that the fractal dimension of the asymp- 
totic cluster (D = 1 in this case) can be extracted from 
growth events infinitesimally close to the initial unit cir- 
cle by decreasing Ao (even without increasing n in this 
specific case). Secondly we consider the deterministic 
itinerary (||) with Golden Mean winding number. Us- 
ing values of F^IO" 6 ) and ^"'(lO^ 5 ) between 1.10 
and 1.20 we can bound the dimension of the cluster to 
1.8305 < D < 1.8380. Note that in this case we need to 
have at least one layer covering which is obtained only 
when n c s» 27r/vAo- Finally, we consider DLA. Here 
the itineraries are stochastic and one c ould imagine that 
only under extensive ensemble averaging one would ob- 
tain tight bounds on D. In fact we find that using values 
of F^IO" 8 ) and F[ n) {2 x 10~ 8 ) between 1.002 and 1.01 
we can bound D as tightly as 1.6896 < D < 1.7135. Note 
that to achieve this accuracy we did not need to go to 

(n) 

high values of F± , but rather used small values of Ao to 
reach convergence very early. This demonstrates again 
the unexpected fact that the asymptotic dimension ap- 
pears as a renormalization exponent right after one or 
a few layers of particles cover the circle, and very much 
before F[ n) ~ n 1 ' . 

At this point we need to address a few questions: 

HI of the 



(n), 



(i) Why classical numerical estimates |11 



fractal dimension of DLA converge so slowly? 

In standard numerical experiments the radius of gyra- 
tion of the grown cluster was plotted in log-log coordi- 
nates against the number of particles, with D estimated 
from the slope. Examining our fixed point functions F£ 
(see Figs. 3-5) we note the slow crossover to linear be- 
haviour, which is not fully achieved even for extremely 
high values of n. In this respect we understand from 
Eq.(^) that a reliable estimate of D from radius of gy- 
ration calculation requires inhuman effort, as was indeed 
experienced by workers in the field ]l2] |. In the present 
formulation the appearance of the asymptotic D as a 



renormalization exponent already at early stages of the 
growth allows a convergent calculation. 

(ii) Is a typical DLA cluster self-averaging? 

It was demonstrated in [0 that an ensemble of DLA 

(n) 

clusters exhibits statistics of F x with standard devia- 
tion that shrinks to zero when n — > oo. Nevertheless in 
p4j it was argued that there may remain residual fluc- 
tuations of the dimension D as extracted from (^). The 
numerics presented above is not sufficient to resolve this 
question, but additional highly accurate numerics in the 
limit Ao — > 0, n — > oo should do. 

(iii) Is the problem co-dimension 1? 

The multi-scaling properties of the harmonic measure 
have left the impression that computing the fractal di- 
mension of DLA will require a simultaneous control of 
the host of exponents characterizing the measure. For 
example the scaling relation D3 = D/2 (with D3 beeing 
a generalized dimension in the Hentschel-Procaccia sense 
pl[) that was derived first by Halsey |2j| strengthened 
this impression. The approach presented here indicates 
that an appropriate fixed point structure can be obtained 
with only one relevant exponent, i.e. 1/D, and this expo- 
nent appears in the dynamics much before the measure 
becomes multiscaling. 
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